Phase Separation Driven by External Fluctuations 
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The influence of external fluctuations in phase separation processes is analysed. These fluctuations 
arise from random variations of an external control parameter. A linear stability analysis of the 
homogeneous state shows that phase separation dynamics can be induced by external noise. The 
spatial structure of the noise is found to have a relevant role in this phenomenon. Numerical 
simulations confirm these results. A comparison with order- disorder noise induced phase transitions 
is also made. 
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The role of noise as an ordering agent has been broadly 
studied in recent years in the context of both tempo- 
ral and spatiotemporal dynamics. In the temporal case, 
which was the one to be addressed earlier, external fluc- 
tuations were found to produce and control transitions 
(known as noise-induced transitions) from monostable 
to bistable stationary distributions in a large variety of 
physical, chemical and biological systems |E|. Spatiotem- 
poral systems have been faced much more recently. In 
these cases, the combined effects of the spatial coupling 
and the noise terms acting upon the system variables may 
produce an ergodicity breaking of a bistable state, lead- 
ing to phase transitions between spatially homogeneous 
and inhomogeneous phases. Results obtained in this field 
include critical-point shifts in standard models of phase 
transitions J^-Q], pure noise-induced phase transitions 
stabilization of propagating fronts §, and noise- 
driven structures in pattern- formation processes [p|-[l2||. 
In all these cases, the qualitative (and somewhat coun- 
terintuitive) effect of noise is to enlarge the domain of ex- 
istence of the ordered phase in parameter space. It is the 
purpose of this Letter to analyse the role of fluctuations 
in a radically different type of spatiotemporal process, 
namely in phase separation dynamics. It will be shown 
that external noise can induce the phase separation pro- 
cess and that this effect is determined by the spatial cor- 
relation of the noise terms. An important conclusion is 



that phase separation does not occur, necessarily, for the 
same range of parameters for which the system presents 
a phase ordering process. 

The dynamics of a large class of spatially-extended sys- 
tems can be described in a general way by the following 
standard model [fili: 



dif)(x, t) 
dt~ 



— if 



£{x,t) 



(1) 



where t) 1S a dynamical field that describes the state 
of the system, J 7 is a free energy functional, and £(x,t) 
is a space-dependent stochastic process that accounts 
for thermal fluctuations. Thermal equilibrium at tem- 
perature T is reached at long times if the fluctuation- 
dissipation relation holds: 



(£(£ , t)£(x', t')) = 2eC(V)5{t - t')6(x - x') 



(2) 



with e = ksT. The system is usually characterised 
by the behavior of the spatial average of the field ip, 
which plays the role of an order parameter. The op- 
erator £(V) does not alter the equilibrium state of the 
system (~ cxp(— J-'/krjT)), but only its transient dynam- 
ics. Two forms of C are usually adopted: for C = 1 (the 
so-called model A) the system evolves towards its equi- 
librium state without any constraint on the value of the 
order parameter; for >C(V) = — V 2 (model B) the order 
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parameter is conserved throughout the dynamical evolu- 
tion. Model A is a prototype of ferromagnetic ordering, 
and model B of phase separation dynamics following a 
quench from a high temperature homogeneous phase to 
a low temperature state. In this latter case, according 
to the value of the order parameter and the quench lo- 
cation, the evolution might proceed either by spinodal 
decomposition or by nucleation |l~4| ]. 

Far from the critical point, the dynamics is not qualita- 
tively affected by the presence of the thermal noise term 
£(x,t). Instead, we will consider the situation in which 
there is an additional source of noise, T](x,t). This hap- 
pens, for instance, when one of the externally controlled 
system parameters is subjected to fluctuations. As these 
fluctuations are external, they do not usually verify the 
fluctuation-dissipation relation (|^) and the system is no 
longer at equilibrium. We will show that, in the case of 
model B, these external fluctuations can induce a phase 
separation in the system. Although this is a reminis- 
cence of noise-induced phase transitions reported earlier 
for model A [pp|JT5[| (and, in what follows, we will make 
a comparison between the effect of external noise in both 
models), an unexpected and notorious feature is that the 
nature of the destabilizing terms due to the noise is in- 
trinsically different for model A and model B. We con- 
clude that phase separation does not necessarily occur 
in the conserved order-parameter model B at the same 
values for which model A shows a noise-induced phase 
transition. It follows that, in contrast with what happens 
at equilibrium, both models have now different station- 
ary distributions. Furthermore, a relevant result of our 
analysis is that the spatial structure of the noise plays 
an important and distinct role for both models. We as- 
sume the following correlation for the external noise with 
a characteristic correlation length A: 

(r 1 (x,t)r 1 (x , ,t'))=2a 2 S(t-t , )g(\x-x , \/X) (3) 

where g is a (short-ranged) spatial correlation function. 
It is expected, and will be confirmed in what follows, that 
since model B represents a domain-growth process, this 
correlation length will play a role far beyond the intuitive 
intensity-reduction effect of space-time noise correlation 
H . Correlation time of the noise should not have a paral- 
lel influence, since all time scales of the system are larger 
than those of the noise. Therefore, it seems that a fi- 
nite (non-zero) correlation time of the noise would not 
differentiate between models A and B, and will not be 
considered here (the time dependence of the correlation 
of r)(x, t) will be assumed to be a Dirac delta, as shown 
in Eq.(|)). 

Although our results are quite general, we shall work, 
for the sake of clarity, with the well known Ginzburg- 
Landau free energy 



T = / dx 



^ 2 + ^ 4 + f|v^| 2 



(4) 



In the absence of additive (thermal) noise, phase separa- 
tion occurs for r < r c — 0. If additive noise is present, 
the transition occurs at r c < 0. We will consider r > 0, 
so that order will not appear spontaneously. The control 
parameter r will be assumed to be subjected to external 
fluctuations, i.e. r — > r + rj(x,t). The spatial correlation 
function g is chosen to be a Gaussian of width A: 
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(5) 



[d is the a?-space dimension) which becomes a delta func- 
tion in the limit A — > 0. It is simpler to analyze the 
role of noise by using a lattice discretization in which the 
space vectors x take values Xi (i — 1, . . . , N), defined on 
regular lattice of linear cell size Ax = 1. The field ip(x, t) 
then becomes a discrete set of variables ipi(t) and simi- 
lar notation is used for the random fields rji(t) and £j(t). 
Under these considerations, the lattice version of model 
(|l|) with the Ginzburg-Landau free energy (jij) is: 



V>; = -C L (ripi + m^i + ^ - KV 2 L ipi) + & 



(6) 



where Cl = 1 for model A and Cl = —^1 f° r model B. 
V| is the lattice Laplacian operator. Finally, in this ver- 
sion, the external noise has a correlation function gu-j\ 
which is the discrete inverse Fourier transform of tjk, the 
corresponding lattice version of the Fourier transform of 
(^), namely (in two dimensions): 

h = cxp (~ (sin(fc x /2) 2 + sin(fc y /2) 2 )) (7) 

There is no closed analytical form for gu—j\ and the de- 
sired values go and g\ (see later) must be obtained nu- 
merically. 

The transition towards an ordered state can be ana- 
lyzed by studying the stability of the homogeneous phase 
ipi = 0. The early time evolution of the statistical mo- 
ments of tpi in Fourier space can be obtained in a linear 
approximation. For example, the second moment (struc- 
ture function) is defined as Sfc(t) = (ipkip-k), where 
N is the number of points of the system. Making use 
of the Stratonovich calculus and Novikov's theorem, its 
evolution equation is [O : 



dSkjt) 
dt 



= -2w(fc)£ fc (t) + ^/(fc) Y,9kS k {t) + 2e 



(8) 



Hence the second moment equation contains a term 
which globally couples Fourier modes and a constant 
term due to thermal noise. The particular values of the 
dispersion relation uj(k) and of the mode-coupling coeffi- 
cient f(k) differ for models A and B. For model A, the 
result is well known Hi 



'eff 



Kk 



f A (k) = 1 



(9) 
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with an effective control parameter r 



>del B the situation is drastically different 



a 2 g . For 



u B {k) = r? ff k 2 



f B (k) = k 2 



(10) 



a 2 V 2 L g and 



with effective control parameter = 
effective diffusion coefficient Kftf = K — o 2 g\. Two main 
differences are observed with respect to model A: the dif- 
fusion coefficient K is also renormalized by the correlated 
external noise, and the noise-induced shift of the control 
parameter r depends now, through the Laplace operator 
V|, on the spatial structure of the noise correlation, i.e. 
not only on the same-site correlation go , but also on the 
nearest-neighbor correlation g\. These differences will 
reveal itself in the position of the transition point where 
the homogeneous state loses stability and phase separa- 
tion appears. 

When neglecting the mode-coupling terms in Eq.(H), it 
is readily seen that perturbations grow when w(k) A '< 
for some interval of k values. We have checked, by 
means of a numerical integration of Eq. (||), that mode- 
coupling terms hardly influence the position of the tran- 
sition curves. Hence the transition point is characterised 
by r e ff — for both model A and model B. The critical 
value of the control parameter r and its dependence on 
the spatial structure of the noise is, however, different in 
the two cases: 



of approximation jjj . In order to corroborate the results 
obtained by means of the linear stability analysis, equa- 
tions (||) have been integrated numerically for models A 
and B in dimension d — 2. A standard stochastic algo- 
rithm is used in order to handle both the additive and 
multiplicative noise terms |16|. Gaussianly-distributcd 
random vectors are generated by means of a numerical 
inversion method, optimised to efficiently produce large 
quantities of Gaussian random numbers [jl7| . According 
to the previous discussion, the spatially correlated exter- 
nal noise is generated in Fourier space with the desired 
correlation function (^), and transformed back to real 



space at each integration time step 11 




model A : 
model B : 



r c 



-<J 2 V 2 L g = o 2 2d{g - 9l ) (11) 



Figure 1 shows the transition curves between homoge- 
nous and inhomogeneous states in the (A, a 2 ) plane for 
models A and B for a fixed value, r = 0.2, of the con- 
trol parameter and spatial dimension d = 2. All points 
located above the curves shown in this phase diagram 
are in an inhomogeneous state, which corresponds to an 
ordered phase in model A (solid curve) and to phase sepa- 
ration in model B (dashed curve). The A— dependence of 
the model- A curve is merely due to the natural "soften- 
ing" effect of noise correlation j|[ll) . In the case of model 
B, on the other hand, additional, non-trivial dependence 
on the correlation length is introduced via the Laplace 
operator. As a consequence, for small A, the transition 
in model B occurs sooner in model B than in model A, 
whereas the situation is the opposite for large values of 
A. A crossing of the two transition curves occurs for an 
intermediate value of A « 1.8. We stress again that the 
presence of ordered regions in the phase diagram of figure 
1 is due to the presence of a multiplicative noise on the 
model, since we are taking r = 0.2 which is larger than 
the mean-field critical value r c = 0. 

The lines drawn in the phase diagram of Fig. 1 have 
been obtained in a linear approximation (11). It is pre- 
sumable that this linear stability analysis will provide 
the position of the transition points up to leading order 



FIG. 1. Phase diagram in the (A, a 2 ) plane for model A 
(solid line) and model B (dashed line), obtained from a linear 
stability analysis (8). Points correspond to numerical simu- 
lations of the full model (6), (diamonds: model A, squares: 
model B) (r = 0.2, K = 1, and e = 1CT 4 ). 

A unified criterion for the existence of an ordered phase 
in model A and of phase separation in model B is the 
growth of the averaged second moment of ip in real space 
(the averaged first moment is not useful for model B be- 
cause this model conserves the order parameter). We 
define this quantity as J(t) = -w (t)) or, alterna- 

tively, as J(t) = Sh(t). The instability point is thus 
defined so that, below it, J(t) decays to a thermal-noise 
background at large times and, above it, it grows to a 
non-zero steady-state value J st . In this way, one can de- 
termine numerically the phase diagram of the system. 
The numerical results are represented in Fig. 1 as dia- 
monds (model A) and squares (model B). It can be ob- 
served that the simulations of the full nonlinear models 
reasonably adjust to the predictions of the linear analy- 
sis. The agreement starts to fail at high values of A and 
a 2 . In fact, according to previous observations in model 
A jl5| and other models , the transition curves might 
be expected to exhibit, for higher values of the noise in- 
tensity, a reentrant branch towards the A = axis. These 
reentrant transitions are hard to observe numerically be- 
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cause of the large noise intensities involved. 

Fig. 2 shows two patterns of a system evolving ac- 
cording to model B, for point 1 in the phase diagram 
of Fig. 1. Depending on the initial conditions we get 
spinodal decomposition (Fig. 2a, with (ip(x, 0)) = 0), or 
nucleation (Fig. 2b, with (i/j(x, 0)) = 0.1). For the same 
values of the noise parameters a homogeneous phase is 
obtained for model A. 
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(b) 

FIG. 2. Spatial patterns of model B: a) Spinodal decom- 
position and b) Nucleation. (t = 2500, A = 0, a 2 — 0.1 and 
e = 10" 4 ) 

For larger values of the noise parameters a reverse sit- 
uation is found. Fig. 3 shows a spatial pattern of model 
A for point 2 in the phase diagram of Fig. 1. Now the 
homogeneous phase corresponds to model B. 




FIG. 3. Spatial pattern of model A. ( t = 100, A = 2, 
a 2 = 6 and e = 10~ 4 .) 

In conclusion, we have demonstrated for the first time 
the ordering role of external noise in processes of phase 
separation. The study has concentrated on the conserved 
time-dependent Ginzburg-Landau model, although our 
results are not restricted to this particular model. Exter- 
nal noise is found to enhance the phase separation pro- 
cess, and this effect is observed to be modified by spatial 
correlation of the noise, which increases the efficiency of 
fluctuations for small correlation lengths and decreases it 
for large correlation lengths. Future work on this issue 
should analyse the important effect of external noise on 
the dynamical scaling of the phase separation process. 
Additional theoretical analysis such as a mean-field ap- 
proach could be useful in determining whether a 
reentrant transition should be expected in this model. 
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